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Abstract 

We observe significant effects of structural distortion on electronic and optical 
properties of CdGa 2 X A (X = S, Se, Te) defect chalcopyrite. The calculation 
is carried out within Density functional theory based tight binding linear 
muffin tin orbital (TB-LMTO) method. Structural parameters and band 
gap of CdGa 2 X 4 agree well with the available experimental values within 
LDA limit. Change in band gap due to structural distortion is 3.63%, 4.0%, 
and 8.8% respectively. We observe significant change in optical properties 
also due to this effect. Effects on optical properties come mainly from optical 
matrix elements. Our results of these response functions agree well with the 
available experimental values. 
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1. Introduction 


CdGa 2 X± (S, Se, Te) defect chalcopyrites are attractive materials because 


, 2[. These compounds 


of their luminescent properties in the visible region 
belong to group II — III 2 — VI 4 defect chalcopyrites. Defect structure is 
due to vacancies at the cation sites. Vacancies are such that they do not 
break translational symmetry of the crystal. Defect structure makes such 
compounds porous. are not only important for device applications 

but they have attracted special attention of the physics community due their 
porisity nature. These vacancies can be occupied by suiatble elements to 
modify their physical properites without breaking symmetry. Therefore it is 
vital to know detailed structural, electronic and optical properties of these 


materials. Severa 


compounds 


expenmen 




al and t 


reoretical studies are reported on these 


id, y y, y y y y y, y . Rang e t. ai. 


9( studied band structure using first principle calculation. But no detailed 


electronic and op 


Zhao et. al. 


ical properties were studied. Recently Ma et. al. 


141 ] and 


191 ] have reported structural, electronic and optical proper¬ 


ties of CdGa 2 Xi (X = S, Se, Te) using DFT based VASP code. It is well 


known that ternary c 


lies 


20, l2l|, l22, |23|, 


ralcopyrites have some interesting structural anoma- 


24] relative to their binary analogous. These structural 


anomalies have significant effects on band gap and so on overall electronic 
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and optical properties. It is well known that p-d hybridization does not have 


significant effect on band gap reduction in the group of II — III 2 — VI 4 chal- 


copyrites 


21 


23], It is a prominant mechanism for band gap reduction in 


other groups of chalcopyrite semiconductors. Therefore the effect of struc¬ 
tural distortion plays a major role in band gap increment/decrement in the 
group of II — III 2 — VI 4 chalcopyrites. But there is no such study, either 
qualitative or quantitative, on these materials. 

Therefore the main motivation in this work is to study effects of structural 
distortion on electronic and optical properties of CdGa 2 X i (X = S, Se, Te). 
We consider here two types of structural distortions, the displacement of an¬ 
ions from their ideal position (u) and structural deformation expressed by 
the ratio of the lattice parameters (c/2a). It has been established that struc¬ 
tural distortion has prominent effects on band gap and electronic properties 


of pure chalcopyrite 


m 


22] and defect/substituted chalcopyrite semiconductors 


24], In the present work we calculate quantitative change in band gap, 


electronic and optical properties due to structural distortion. Structural pa¬ 
rameters like lattice parameters, tetragonal distortion, anion displacement 
parameters and bond lengths are calculated by structural relaxation using 
energy minimization proceedure. There have been no calculation so far of 
bulk modulus of CdGa 2 Te 4 . We calculate bulk modulus by fitting our data 
of energy versus volume to Murnaghan equation of states 25]. We carry 
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out detailed electronic properties, like band structure, total density of states 
(TDOS) and partial density of states (PDOS) and then calculate effects of 
structural distortion on TDOS. To study effects of structural distortion on 
optical properties we calculate joint density of states (JDOS) and square of 
optical matrix elemnets (OME) separately. With these calculations we show 


how structural distortion have effects on optical response fune 
dilations are based on density functional theory (DFT) 


26 


ions. Our cal- 


27] first principle 


tec 


inique, tight binding linearized muffin-tin orbital (TB-LMTO) method 


28 


29| . Basis functions in TB-LMTO method are localized. Therefore, 


very few such functions are required to represent highly localized d-orbitals. 
Hence the calculation is not only cost effective, it gives also accurate results. 
DFT requires a suitable choice of exchange and correlation functional for 
the electron-electron correlations. We have chosen such a functional within 


LDA(Local density approximation) 30j. It is known that LDA underesti¬ 
mates band gap in semiconductors. Therefore it is necessary to use more 


advanced correlation functional like hybrid functional 


311 ]. But our moti¬ 


vation here is not to calculate band gap accurately. Purpose is to show 
effects of structural distortion on electronic and optical properties. It is 


bund 


22 


23 


3 A gives good structural properties of chalcopyrite semiconductors 


24 


3,2] and calculate accurate band gap reduction due to p-d hy¬ 


bridization. This is because reduction does not depend significantly on the 
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kind of functional used [32]. We expect the same to be true also for reduction 
in band gap due to structural distortion. 


2. Methodology 

2.1. Structural & Electronic Properties 

The ab initio method is based on Density Functional Theory of Kohn- 
Sham 


26 


27]. The total electronic energy is a function of electron density 


which is calculated using variational principle. This requires self-consistent 
calculations. In practice the Kohn-Sham orbitals are usually expanded in 
terms of some chosen basis function. We have used the well-established TB- 


LMTO-ASA method, discussed in detail elsewhere 


28 


for the choice of 


the basis function. To solve the one-electron Schrodinger-like equation self- 
consistently, we use the scalar-relativistic linear muffin-tin orbital method in 
the atomic sphere approximation (ASA) which includes the combined correc¬ 
tion terms. The Hamiltonian and overlap matrices are evaluated over atomic 
spheres where the potential is spherically symmetric. Electron correlations 


29j ]. We have used the von BarthHedin ex- 


are taken within LDA of DFT 
change 30] with 512 k-points in the irreducible part of the Brillouin zone. 
The basis of the TB-LMTO starts from the minimal set of muffin-tin orbitals 
of a KKR formalism and then linearizes it by expanding around a ’nodal’ 
energy point represented by EKAP in band structure figures. The wave- 
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function is then expanded in this basis. Basis functions np to l ma x= 3 for 
Cd and l m ax= 2 for both Ga & X (S, Se, Te) are used for the band structure 
calculation. Cd 4d, 5s & 5p, Ga 4s & 4p, S 3s & 3p, Se 4s & 4p, Te 5s & 
5p and empty spheres (E, El, E2, E3) Is linearized muffin tin orbitals are 
treated as valence states. The Cd 4f Ga 4d, S 3d, Se 4d, Te 5d & 4f and 
empty spheres p and d states are downfolded for all the three systems in our 
calculation. These downfolded states are included in the tails of the valence 
LMTO’s. These states do not contribute to the dimension of the Hamiltonian 
(H) and the overlap (O) matrices, but carry charge. This technique avoides 
the problem of ’ghost’ bands. In our calculaion the core states are treated 
relativistically, while the valence states were calculated scalar relativistically. 
Defect chalcopyrite semiconductors are loosely packed crystals (low packing 
fraction). To overcome such unfavorable physical situation, it is necessary 
to introduce empty spheres at appropriate chosen interstitial sites for self- 
consistent calculation. This technique is incorporated in the advanced TB- 
LMTO-ASA method. We ensure proper overlap of muffin tin spheres for 
self-consistency and the percentage of overlaps are found to be 8.64%, 8.41% 
and 8.40% for CdGa 2 X 4 (X = S, Se, Te) respectively. The sphere radii and 
the positions are chosen in such a way that space filling in the unit cell is 
achieved without exceeding a sphere overlap of 16%. The atomic and inter¬ 
stitial sphere radii are calculated and given in table 1, table 2 and table 3 
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respectively for CdGa 2 X 4 . Total nine empty spheres of four different radii 
(E, El, E2 and E3) are introduced for proper overlap of the spheres in all 
three systems. We obtain the ground state equilibrium values of lattice pa¬ 
rameters, anion positions (u) and bond lengths by structural relaxation. It 
is done by minimizing total energy using the above first principle procedure. 
We calculate the Bulk modulus by hiting energies versus cell volume data to 
Murnaghan equation of state. 

2.2. Optical Properties 

Optical properties can be determined by the dielectric function e(u ;) which 
is a complex quantity. e(uj) = £i(co) + i£ 2 (uj) where £i(cu) & £ 2 (co) are the real 
and imaginary parts of the dielectric function. There are two contribution 
to £(oj). One is due to intraband transitions and another is due to interband 
transitions. The contribution from intraband transitions is important for 
metals and semimetals. But interband transitions are important in case 
of semiconductors. To calculate the direct interband contributions to the 
imaginary part of the dielectric function £ 2 (uj), all the possible transitions 
from the occupied to the unoccupied states within the irreducible Brillouin 
zone must be summed up. The interband transitions can be of two types 
i.e. direct and indirect transitions. Here we neglect the indirect interband 
transitions which involve electron phonon scattering and are expected to give 
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Table 1: Radii and positions of atomic and empty spheres in CdGa^Sit. 


atoms 

X 

y 

z 

sphere radii (a.u.) 

Cd 

0.000 

0.000 

0.000 

2.828 

E 

0.000 

0.500 

-0.497 

2.332 

Gal 

0.000 

0.000 

0.995 

2.585 

Ga2 

0.500 

0.000 

-0.497 

2.642 

S 

0.248 

-0.270 

-0.265 

2.572 

S 

-0.248 

0.270 

-0.265 

2.572 

s 

0.230 

0.252 

-0.730 

2.572 

s 

-0.230 

-0.252 

1.259 

2.572 

El 

0.500 

0.000 

-0.041 

2.582 

El 

0.500 

0.000 

-0.953 

2.582 

E2 

0.270 

-0.248 

0.265 

2.572 

E2 

-0.270 

0.248 

0.265 

2.572 

E2 

-0.252 

-0.230 

0.730 

2.572 

E2 

0.252 

0.230 

-1.259 

2.572 

E3 

0.000 

0.000 

0.456 

2.367 

E3 

0.000 

0.000 

-0.456 

2.367 





Table 2: Radii and positions of atomic and empty spheres in CdGa, 2 Se 4. 


atoms 

X 

y 

z 

sphere radii (a. 11 .) 

Cd 

0.000 

0.000 

0.000 

2.827 

E 

0.000 

0.500 

-0.499 

2.400 

Gal 

0.000 

0.000 

0.997 

2.646 

Ga2 

0.500 

0.000 

-0.499 

2.618 

Se 

0.249 

-0.264 

-0.265 

2.727 

Se 

-0.249 

0.264 

-0.265 

2.727 

Se 

0.236 

0.251 

-0.732 

2.727 

Se 

-0.236 

-0.251 

1.263 

2.727 

El 

0.264 

-0.249 

0.265 

2.726 

El 

-0.264 

0.249 

0.265 

2.726 

El 

-0.251 

-0.236 

0.732 

2.726 

El 

0.251 

0.236 

-1.263 

2.726 

E2 

0.500 

0.000 

-0.028 

2.647 

E2 

0.500 

0.000 

-0.970 

2.647 

E3 

0.000 

0.000 

0.471 

2.427 

E3 

0.000 

0.000 

-0.471 

2.427 
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Table 3: Radii and positions of atomic and empty spheres in 


atoms 

X 

y 

z 

sphere radii (a.u.) 

Cd 

0.000 

0.000 

0.000 

2.895 

E 

0.000 

0.500 

-0.494 

2.589 

Gal 

0.000 

0.000 

0.989 

2.666 

Ga2 

0.500 

0.000 

-0.494 

2.791 

Te 

0.251 

-0.264 

-0.253 

3.016 

Te 

-0.251 

0.264 

-0.253 

3.016 

Te 

0.236 

0.249 

-0.736 

3.016 

Te 

-0.236 

-0.249 

1.242 

3.016 

El 

0.264 

-0.251 

0.253 

3.015 

El 

-0.264 

0.251 

0.253 

3.015 

El 

-0.249 

-0.236 

0.736 

3.015 

El 

0.249 

0.236 

-1.242 

3.015 

E2 

0.500 

0.000 

0.000 

2.681 

E2 

0.000 

0.500 

0.000 

2.681 

E3 

0.000 

0.000 

0.467 

2.613 

E3 

0.000 

0.000 

-0.467 

2.613 
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only a small contribution to s(co). The expression for the imaginary part of 
the dielectric function is g iven by the well known Kubo formula for the direct 


interband transitions 


33, 


34] 


e]M = = 


—87T e 1 
2m* 2 fl u 2 


E E ' Pl*)|Vi(l - F,)S(E, -E t - hw) (1) 


i / 


where, m* is the effective mass of the electron,£7 is the volume of the sample, 
i and f refer to the initial and final states respectively, 7 refers to the direction 
of polarization of the incoming photon, p is the momentum of the electron and 
Fi f is the occupation probability of the initial and final states respectively. 
For semiconductors, i lies in the valence band and f in the conduction band 
and at OK we have Fi = 1 and Ff = 0. We can obtain the real part of the 


dielectric function ei(cu) from Kramers-Kronig relation 


33]. 


/ . 2 r 00 (co' - co)e 2 (co') , 

ei(o;)-1 + H-/ - -2 -5 —du 

tt Jo u — co z 


( 2 ) 


A number of methods have been proposed 


within the frame work of the LMTO 
and semiconductors 


Hobbs et.al. 


35, 


3(: 


or calculating optical properties 


m 


38] for both metals 


35 


37]. I 11 this study we use the method developed by 


which avoids the determination of gradient operator. This 


method allows for the inclusion of non-local potentials also in the Hamil¬ 
tonian. The momentum matrices which appear in the Kubo formula, are 
written in terms of Gaunt coefficient and potential parameters which are de¬ 


fined within the TB-LMTO method. We calculate these Gaunt coefficient 
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separately. 


The imaginary part of dielectric function, optical matrix elements and the 
joint density of states can be calculated by performing an integration over 
the Brillouin zone (k-space integration). For all these calculation, the tetra¬ 
hedron method is used within TB-LMTO formalism. In this method, the di¬ 
electric function is expressed as an integral over the constant-energy surface 
, E c (k ) — E v (k). The eigen values and eigen vectors are then calculated on a 
mesh in the irreducible Brillouin zone. This zone is devided in to tetrahedra 
of equal volume and the mesh of k-points (maximum 512 points) defines the 
corners of each tetrahedron. The interpolated function is continious at the 
boundaries of the tetrahedra. The irreducible Brillouin zone is completely 
divided into tetrahedra. To calculate different necessary formulae for optical 
properties, the output data from the TB-LMTO band structure calculations 
are used. 


3. Result and discussion 

3.1. Structural properties : 

One unit cell of CdGa 2 X 4 chalcopyrite is shown in figure 1. It is basicaaly 
a supercell of two zinc-bende unit cells placed one on top of the other with 
a 90° rotation about c-axis. There are two Cd, two vacancies, four Ga and 
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eight X atoms per unit cell. The positions of various atoms in the tetragonal 
unit cell are (in Wyckoff notation): Cd on 2a site (0 0 0), Gal on 2b site (0 0 
0.5), Ga2 on 2d site (0 0.5 0.75), vacancy on 2c site (0 0.5 0.25) and X on 8g 
site [u x u y u z ). Here u x , u y and u z are anion displacement parameters along 
three axes. The space group of the defect chalcopyrite is IA(Sf). Unit cell 
with rj{= c/2a) = 1, u x , u y and u z equal to 0.25, 0.25 and 0T25 respectively, 


is referred as ideal structure (ideal case) 


20 


, 2jJ, 22, [231 24], Anion shifts 


along all the three directions in these systems, unlike only along x-direction 


found in the case of pure chalcopyrites 


20 


22], This new position of anion 


is called anion displacement. This is due to the reduction in symmetry in 
these systems. Therefore, all cations-anion bond lengths are inequivalent in 
these systems. Calculated bond-legnths tabulated in the table 4 show an¬ 
ion acquires an equilibrium position closer to the vacancy than to the other 
cations as expected. 

Table 5 shows our result of lattice parameters a & c, tetragonal distortion, 
and anion displacement agree quite well with experimental and other calcu¬ 
lated results. Table 6 shows our result of bulk moduli agree well with the 


other theoretical and experimental results also. 
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Table 4: Bond lengths(R) in A 


CdGa 2 Si 

Present Work 

Experimental 0 

Other Theoretical 6 

Rcd-s 

2.48 

2.52 

2.57 

R-Cal-S 

2.37 

2.29 

2.31 

RGa2-S 

2.40 

2.33 

2.31 

Ry acancy—S 

2.26 

- 

- 

CdGa 2 Se4 

Present Work 

Experimental 0 

Other Theoretical b 

Rcd-Se 

2.57 

2.63 

2.70 

Rcal-Se 

2.48 

2.41 

2.41 

RGa2-Se 

2.47 

2.41 

2.41 

Ry acancy—Se 

2.37 

- 

- 

CdGa 2 Te 4 

Present Work 

Experimental^ 

Other Theoretical 6 

Rcd-x 

2.74 

- 

2.88 

Rcai-x 

2.63 

- 

2.68 

Rca2-X 

2.69 

- 

2.66 

Ry acancy—X 

2.59 

- 

- 


a Ref. [10]; b Ref. [14]; c Ref. [39]; d Ref. [40]; e Ref. [19] 
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Figure 1: One unit cell of CdGa 2 X 4 chalcopyrie semiconductor (a) front view (b) top 
view, (c) primitive cell of CdGa 2 Te 4 . 


Table 5: Structural parameters. 


Compounds 


a 

(A) 

c 

(A) 

^X 

u y 

u z 

CdGa 2 S A 

Present work 

5.51 

10.94 

0.270 

0.248 

0.133 


Experiment® 

5.53 

10.16 

0.271 

0.261 

0.140 


Other theory 6 

5.64 

10.34 

0.268 

0.271 

0.136 

CdGa 2 Se4 

Present work 

5.72 

11.41 

0.264 

0.249 

0.133 


Experiment® 

5.73 

10.70 

0.275 

259 

139 


Other theory 6 

5.88 

10.78 

0.273 

0.270 

0.138 

CdGa 2 Te4 

Present work 

6.18 

12.22 

0.264 

0.251 

0.128 


Experiment^ 

6.10 

11.70 

0.270 

0.260 

0.135 


Other theory 6 

6.27 

12.00 

0.272 

.255 

137 


a Ref. [10]; 6 Ref. [14]; c Ref. [39]; d Ref. [40]; e Ref. [19] 
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Table 6: Bulk modulus B in GPa. 


Systems 

Present work 

Experimental 

Other theoretical 

CdGa 2 Si 

57.10 

64.0 a 

58.44 6 

CdGa 2 Se4 

38.32 

41.0 C 

36.13 6 

CdGa 2 Te± 

31.00 

- 

- 


a Ref. [10]; b Ref. [14]; c Ref. [8] 

3.2. Structural effect on electronic properties 

Band structure, total density of states (TDOS) and partial density of 
states (PDOS) are calculated to study the detailed electronic properties of 
CdGa 2 X 4. The Fermi level is set at zero energy in the band structure and 
total DOS plots and marked by dashed horizontal line. E, El, E2 and E3 
represent empty spheres of different radii in band structure diagram which 
are introduced at appropriately chosen interstitial sites. Band strucrure di¬ 
agram of CdGa 2 S^ for non-ideal structure is only shown here in the figure 2 
for reference. A similar picture arise in the case of other two systems. The 
partial density of states shown in the figure 2 show that p-d hybridization 
play important role in band formation. 

The valence band maximum and conduction band minimum are located 
at the center of brillouin zone denoted as “G” (T point). This indicates that 
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these are direct band gap semiconductors. Table 7 shows our result of band 
gaps underestimate the experimental result and are closed to some of the 
other theoretical results. This die to LDA. It is known that LDA underesti¬ 


mates banc 
like hybrid 


gap by 30-50% 4l|. There are are better functionals avialable 


311 ] for accurate calculation of band gap. But our motivation here 


is not to calculate accuarate band gap. Our main motivation is to show ef¬ 
fects of structural distortion on electronic and optical properties. The choice 
of LDA does not affect such a study. 

We had earlier carried out quantitative study of effect of structural dis¬ 
tortion on band gap in a series of chalcopyrite semiconducting compounds 


22 


23 


24], We apply similar method here in the present study also. Table 8 


shows % of increment in band gap due to the structural distortion. Figures 
2 & 3 show no significant differences in the structures of TDOS in ideal and 
non-ideal cases. A very little differences are visible in the conduction bands 
in the cases of CdGa, 2 S .4 & CdGa 2 Se 4. The conduction band width increases 
in case of ideal compared to the non-ideal case. Whereas signifi¬ 

cant differences are visible in both valence band and conduction bands for 
ideal and non-ideal GdGa 2 Te4. We also observe that the effect of structural 


distortion on band gap increases from CdGa 2 S' 4 to GdGa 2 Te 4 . 
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Figure 2: (Top panel): Band structure and density of states for non-ideal CdGa, 2 Si.. 
(Bottom panel): Density of states for non-ideal & CdGa 2 Te 4 
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Figure 3: Total DOS for ideal CdGa 2 X 4 (X = S,Se,Te ) 


Table 7: Energy band gap E g (eV). 


Systems 

E g 

rp exp a 
^ 9 

rpother 

E 9 


(eV) 

(eV) 

(eV 

CdGa 2 S 4 

2.20 

3.25-3.44 

3.30 fc , 2.11 c 

CdGa 2 Se 4 

1.75 

2.57 

2.76 fe , 1.46 c 

CdGa 2 Te 4 

1.25 

1.50 

1.90 fe 


E e q xp : experimental, E° ther : other calculated gap 


Ref. 


42], b Ref. 


9J; c Ref. 


14]; 
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Table 8: Effect of structural distortion on band gap (eV). 


Systems 

Ideal 

Non-ideal 

Increment in 


(eV) 

(eV) 

band gap (%) 

CdGa 2 S 4 

2.12 

2.20 

3.63 

CdGa 2 Se A 

1.68 

1.75 

4.0 

CdGa 2 Te4 : 

1.14 

1.25 

8.8 


3.3. Structural effects on optical properties 

There are two important factors in the Kubo formula for imaginary part 
of dielectric function e 2 (tu). One is joint density of states (JDOS) and the 
other is the “square of the optical matrix elements” (OME) (transition prob¬ 
abilities or electron photon interaction matrix elements). The various peaks 
in 62 ( 01 ) come from the joint density of states for values of E c & E v at which 
slopes of valence band and conduction band are equal (Von ffove singular¬ 
ities). This occurs at high symmetry points in the Brillounin zone. For 
other values where slope of E c is not equal to slope of E v , transitions may 
occur at any general points in Brillouin zones. But OME is a difficult factor 
to be interpreted accurately. In almost all the works/literature, there is no 
explicit calculation of this quantity to see its effect on optical absorption 
of chalcopyrite semiconductors. The observed or calculated absorption co¬ 
efficient/imaginary part of the dielectric function are generally interpreted 
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as transition from various critical points. Therefore majority of the work 
correlates 62 with joint density of states. 

Chalcopyrite semiconductors have tetragonal structure and therefore its 
optical properities show anisotropic nature. Therefore oprical response func¬ 
tions are different when incoming photon is polarised along c-axis than when 
it is polarised along perpendicular to c-axis. 

Figure 4 shows no significant differences in JDOS between ideal and non¬ 
ideal structures of The magnitude of JDOS is almost equal in both 

the structures . Similar results are also observed in other two systems. But 
figure 5 shows remarkable differences in OME between the two structures. 
The magnitude of optical matrix elements is almost 1.75 times higher in 
ideal case than the non-ideal. Lot of structures which are present in the 
ideal structures get drastically reduced in non-ideal structure. Therefore any 
differences in 62 between the two structures are mainly coming from changes 
in OME. For example we notice from figure 6 that the magnitude of the 
imaginary part of dielectric function ( 62 ) in the case of ideal CdGci 2 S 4 is 
almost double compared to non-ideal The highest peak in OME 

and e 2 is found at the same energy. In case of ideal CdGa 2 Te 4 the magnitude 
of OME increases up to 7 and the peak found at energy 2.5 eV in non-ideal 
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Figure 4: (Left panel) : Joint density of staes 
CdGa 2 S 4 



for non-ideal; (right panel): for ideal 


case splits in to two sharp peaks. So structures in OME are clearly reflected 
in e 2 for both CdGa 2 Se 4 and CdGa 2 Te 4 . Significant differences are also 
observed in the case of real part of dielectric function as can be seen in 
the figure 7. Figures 5, 6 & 7 show also anisotropic nature of the optical 
properties and they are more prominent for the non-ideal structure than the 
ideal. 

We also find that static dielectric constants and refractive indices increase for 
ideal structure of CdGa 2 X 4 . Table 9 & 10 show the static dielectric constant 
€i (0) and refractive index n(0) for non-ideal and ideal CdGa 2 A" 4 respectively. 
Our calculated 6 \ and e 2 for C'dGa 2 Te 4 agree well with the result of Ozaki 
et.al. 45] 
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Figure 5: (Left panel) : Optical matrix elements for non-ideal (right panel) : for ideal 


CdGci 2 X 4 ; X = S,Se,Te respectively 
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Figure 6: (Left panel) : Imaginary part of dielectric function for non-ideal (right panel) : 


for ideal CdGd 2 X± ; X = S,Se,Te respectively 
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Figure 7: (Left panel) : Real part of dielectric function for non-ideal (right panel) : for 


ideal CdGa, 2 X 4 ; X = S,Se,Te respectively 





























Table 9: Static dielectric constant ei(0) and refractive index n(0). for non-ideal 


Optical Constants 

CdGa 2 S 4 

CdGa 2 Se 4 

CdGa 2 Te 4 

0(0)^ 

7.32 

7.61 

10.94 

ei(0)H 

5.35 

5.84 

7.71 

O ( Smother 

4.29“ 

5.56“ 

9.53 6 

e l(0 ) other 

4.22“ 

5.45“ 

9.56 b 

ei(0 )expt 

- 

- 

12.04 c 

O(O)expt 

- 

- 

11.47“ 

tRo)- 1 

2.70 

2.75 

3.30 

77,(0) 11 

2.31 

2.41 

2.77 

Mother 

2.07“ 

2.35“ 

3.09 6 

^(0 ) other 

2.05“ 

2.33“ 

3.09 6 

n (°)ix P t 

2.3 d 

- 

- 


- 

- 

- 


a Ref. [14]; b Ref. [19]; c Ref. [43]; c Ref. [44] 

Table 10: Static dielectric constant eRO) and refractive index n{ 0) for ideal CdGa 2 X 4 . 


System 

o(0)^ 

ei(O) 11 

n(O)- 1 

77,(0) II 

CdGa 2 S 4 

11.09 

12.89 

3.33 

3.59 

CdGa 2 Se 4 

11.65 

15.32 

3.41 

3.91 

CdGa 2 Te 4 

15.95 

14.96 

3.99 

3.86 
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4. Conclusion 


Calculations and study of structural, electronic and optical properties of 
CdGa 2 X 4 (X = S,Se,Te ) suggest that all these three compounds are direct 
band gap semiconductors with band gap 2.20, 1.75 and 1.25 eV respectively. 
Our study further shows that electronic properties such as band structure, 
TDOS and PDOS of these semiconductors significantly depend on the struc¬ 
tural distortion. All these calculations are carried out using DFT based 
TB-LMTO method. We use LDA for our exchange corelation functional. 
Taking into account the underestimation of band gap by LDA, our result 
of band gap and structural properties agree with the available experimental 
values. We find an inecreament of band gaps by 3.63%, 4.00% and 8.8% 
respectively for CdGa 2 Xi (S, Se, Te) due to structural distortion. We also 
find significant effects of structural distortion on optical properties of these 
compounds. We have calculated joint density of states (JDOS), square of 
optical matrix elements (OME) and the optical response functions like real 
and imaginary part of dielectric functions, static dielectric constants. We 
conclude that main effects of structural distortion on optical properties come 
from the significant changes in OME due to structural variations. Results of 
real and imaginary parts of dielectric constants, static dielectric constants of 
these compounds agree well with the available experimental results. 
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